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Abstract 



Completeness of the eigenfunctions of a quantum mechanical system is crucial for the 
probability interpretation. By using the method of contour integral we give proper nor- 
malized eigenfunctions for both discrete and continuum spectrum of the Morse potential, 
and explicitly prove the completeness relation. As an application we use our spectral 
decomposition formula to study the problem of the pricing of an Asian option traded in 
the financial market. 



1 Introduction 



Since its first introduction [Tj to quantum physics in 1929, the morse potential as a 
solvable model has been studied extensively by many different approaches, see e.g. [2]-[H]- 
Most of these works focus only on the discrete spectrum. In this paper we will discuss 
the completeness issue of its eigenf unctions for both discrete and continuum spectrum. 
Completeness is a very important property of a quantum mechanical system. It is crucial 
for the probability interpretation. To establish the completeness of a set of eigenfunctions, 
it is necessary to properly choose their normalization constants, which will give the correct 
integral measure in the spectral decomposition of a square-integrable wave function. The 
main purpose of the present paper is to do so for the Morse potential, especially for its 
continuum eigenstates. 

There is a very powerful method [TO], due to E. C. Titchmarsh, to study the above- 
mentioned spectral expansion problem for a given self-adjoint Sturm-Liouville problem. 
The advantage of this method is the fact that it uses only classical complex analysis, 
which is rather familiar to physicists, and does not require the knowledge of the abstract 
operator theory in functional analysis. By using Cauchy contour integrals, the Titchmarsh 
method will automatically gives a complete system of proper normalized eigenfunctions. 
The poles inside the contour correspond to the discrete spectrum, while the contributions 
from the cuts represent the continuum spectrum. It can even deal with subtleties related 
with self-adjoint extensions. 

In this paper we will apply the Titchmarsh method to establish the completeness 
of the eigenfunctions of the Hamiltonian with Morse potential. In section 2 we quickly 
review some main points of the Titchmarsh expansion theory. In section 3 we apply this 
method to the Morse potential. We give the proper normalization for both discrete and 
continuum eigenfunctions, and prove the completeness relation explicitly. In section 4 we 
discuss an interesting application of our result to the pricing of an Asian option in the 
financial market. 

2 Review of the Titchmarch expansion theory 

In this section we will review some main points of the Titchmarsh expansion theory [10] . 
For a stationary Schrodinger-type equation with a real potential function V(x) 
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we want to get the eigenvalue A and the corresponding eigenfunction ip\(x) satisfying some 
proper boundary conditions. The discrete eigenvalues are interpreted as energy levels of 
possible bound states, while the continuum part corresponds to scattering states. For 
the consistency of the probability interpretation of quantum mechanics, it is crucial that 
the set of all eigenfunctions forms a complete basis of the Hilbert space. Without this 
property, probability is not conserved and unitarity is lost. If x G (a, b) with finite a, b and 
the potential function V(x) is regular at both endpoints, this problem is called regular. 
However if one of or both of two points tend to infinity, or the potential becomes singular 
there, the problem is called singular. Since we want to apply this method to the Morse 
potential, we assume that x G (—00, +00). 

Define functions <j)(x, A) and 8(x, A) to be solutions of (12. ip with boundary conditions 

0(O,A) = O, 0'(O,A) = -1; (2.2) 
0(O,A) = 1, 0'(O,A)=O. (2.3) 

We introduce a finite cutoff b > instead of +00, and consider the following mapping 

, , AX 9(b, X)z + d'Cb, A) , , 

^ ( ' ;M):= - *(mU»-(m) • (2 ' 4) 

This fractional linear transformation maps the real axis Imz = to a circle. It is proved 
by H. Weyl that, under the limit b — >■ +00, the image of real axis in the z-plane either 
contracts to a limit point or converges to a limit circle. The Weyl-Titchmarsh m-function 
m + (A) to value of the limit point or any point on the limit circle. Then it follows that 

ip + (x, A) = 6(x, A) + m+(A) (f)(x, A) (2.5) 

is a solution of (12. ip and normalizable at +00. Similarly, there is another m-function 
m_(A) for the left endpoint such that 

1>-(x, A) = 6(x, A) + m_(A) </>(x, A) (2.6) 

is normalizable at —00. It can also be shown that, if in the limit point case, these are the 
only normalizable solutions for each endpoint. Define the resolvent kernel G(x,x'; A) as 

G{x,x;X) = — — — , [x>x ); 

m + (X) — m_(Aj 



•4>-{x, A) ij} + (x, X) 
m + (A) — m_(A) 



(x < x) . (2.7) 



It is the unique solution of 

d 2 G 



(J)2 +(V(x)-X)G= 8(x-x'), (2.8) 
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Figure 1: The integration contour C := C\ + L + + Cg + L- which is used in (|2.1ip and f|3.20[) . 
We assume that A = is a branch point, and on the negative real axis A < there are some 
poles. We use the convention about the argument: argA = for A G L_. 

which is symmetric between x and x' and square-integrable with respect to each variable. 
For any wave function f(x) to be expanded, let 



The function A) can be analytically continued to complex A with possible poles and 

branch points. For the sake of simplicity and further applications in the next section, we 
assume the function 3>f(x; A) is meromorphic in the complex A-plane cut open along the 
positive real axis A > 0. By Cauchy's theorem we have 




(2.9) 



which is the solution of the inhomogeneous equation 




(2.10) 
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where the contour C is chosen as shown in Figure [TJ For reasonable potential functions 
it can be proved that 

^- (f A) d\ ->■ f(x) , A^oo. (2.12) 
2tt« Jc a 

Therefore we obtain finally an equation as 

1 f°° 

f(x) = -J2 Res$ /+ — J (<f> f (x;\e- 2 ™)-<!> f (x;\))d\. (2.13) 

poles ^ 

This is actually the desired expansion formula for wave function f(x), and the complete- 
ness follows immediately. From this quick review we can see the advantage of this method. 
It does not require any abstract operator theories, and everything can be obtained just 
by the tools of classical complex analysis. So it is especially convenient for physicists. 



3 Spectral decomposition of the Morse potential 

Suppose k > 0, the Morse potential is pQ 

V M (x) = k 2 ( e - 2{x - Xo) - 2 e - (x - Xo) ) . (3.1) 

The corresponding Schrodinger equation is 

d 2 ^ 



dx 2 



+ V M (x)^(x) = Xip(x). (3.2) 



Define 



we have 



u 



2 Ke - {x - xo) , w(u) = Vuij)(x(u)) , (3.3) 



d 2 w (Ik 1/4 - (-A) 1 , , 

This is the Whittaker equation, so we have 

w(u) = Cl M K ^ (u) + c 2 W K ^ x (it) . (3.5) 
Therefore the solution of ip is 

i/j(x) = cii>i(x) + c 2 i>2(x) 

= ^e^^^M^^e-^+c.W (3.6) 



The Wronskian W(ipi,ip2) is 



W{ikM = Mx)ti(x)-&(x)Mx) = y (1 + 2?v ^ . (3.7) 

r( ± — k + iy \ ) 



When x — > — oo, i.e. u — > +00, we have 

* 1 



, = -=M ~ (a + e»/ 2 «-« + a_e-«/V) , (3.8) 

V" 

^2 = 4t^^(«) ~ e "" /V - (3-9) 



When x — > +00, i.e. w — > 0, we have 



«'l = -^M^t*) ~ "' vA - (3-10) 



^ = -^^^(u) ~ + . (3.11) 



If — 2tt < argA < 0, then Im\/A < 0, so the solution normalizable at —00 is ^(x), while 
the solution normalizable at +00 is ip\{x). Therefore both singular points are limit points. 
The resolvent kernel G(x,x'; A) is the normalizable solution of 

- + (^(*) - A) G = 5(x - x') , (3.12) 

ax z 

which is continuous at x = x' but 



dG 

dx 



X=x'+£ 



dG 

dx 



-1, e^0 + . (3.13) 



It can be checked that the solution ia3 



G(x,x';X) = G + (x,x';X) = * + n^/aO ».»v^( ^ x > x' ; 

r(l + 2zVA) V« vV 

= G_(x, x ; A) = = —= —= , x < x . 

T{l + 2iVX) Vti 

This resolvent kernel G(x,x'; A) has simple poles at 

A n = -( K - n - 1/2) 2 , n = 0,l,--- ,[k-1/2). (3.14) 



1 In this paper we choose — 2tt < argA < 0. If we use the usual convention < argA < 2tt, then iy/X 
should be replaced by —i\X everywhere. These two different choices will give us the same normalized 
eigenfunctions as in (|3.28l) and (I3.29[) . 
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which come form the factor T(l/2 — k + iy/X) since^l 

r(l/2 ft -f- ) = 4^ + 0(1) . (3.15) 

n! A — \ n 

Note that the upper limit k — 1/2 of n follows from the requirement Ihia/A < 0, which is 
a direct sequence of — 2tt < argA < 0% At these poles we have 

M K>iVX -Ju) = n r ( ( 2 ^r n 2 ) n) ^I^) , (3.16) 

= (-^"nle^/V^Lf- 2 "- 1 ^), (3.17) 

where is the generalized Laguerre polynomial. Therefore the residues of G(x,x'; A) 
at these poles are 

"! <^ 2 " 1) e - u /2 MK -n-l/2 L 2 K -2n-l (M) g-u'^, «-„-l/2 L 2«-2n-l (3 lg) 



r(2«-n) 

For any function f(x), the solution $/(x; A) of $" + (A — V(x)) $ = /(a;) is 



$/(z;A) = / a/; A) /(a/) dz' 

G+{x,a?;\)f{a?)da? + / G_(x, a/; A) f(x') dx' . (3.19) 



As reviewed in the previous section, to obtain the spectral decomposition of the function 
f(x), we need to consider § c &f(x; A) dX with the contour C defined in Figured! The 
Cauchy theorem tells us 



& $ f (x; A) dX + 2iri V Res $/ = . (3.20) 
Jc , 

poles 

The contributions from poles, as calculated above, are 



2m nl %~ 2 ^ 1] e-^u^-^Ll^-'iu) 



n=0 

X 



/oo 
e - u7 2 w /«- n -l/2 L 2 re -2n-l( u /) /( y) ^ _ (33^ 
-00 



2 Note that the other factor M K i r^(u)/T(l + 2i\X) has no poles in the A-plane. 

3 This upper limit does not depend on the convention about the argument range of A. If we choose 



< A < 2vr, it still holds. 



6 



Next we consider the sum of two integrals along the positive axis 



h+ + h 



L+ 
x 



$f(x;\)d\+ J $/(x;A)cfA 

POO 

dx'f(x') / dX (G+(x,x';X) -G+(x,x';Xe- 2nl )) 
Jo 

■OO POO 

+ / dx'f(x') / dX (G-(x,x';X) -G-(x,x';Xe- !M )) 



JO 

OO 



U 



oo Jo I T(l + 2iy/\) 



r(±-K-iVX)M K> _ iV - x (u)\ W K ^{u>) 



r(l-2i v / A) V* 1 



v! 



u 



Jo 1 r(l + 2^ v / A) \/" 



ra- K -iVx) M Kt _ isrx (u')\ w K>iV - x (u) 



r(i-2iVA) 



1 r 00 



dA sinli(27rVA) 



1 



2 



1 ( ^-k + zv^ 



x 



Here we have used the relation W K -n(u) = W K ^(u) and 



f(x') dx' 



r(-2/x) 



r(2//) 



r (| — « — //) r(| — « + //) 

By (2.15.2) of [10], for square-integrable function f(x), we have 



$ f (x;X) 



+ O 



|A|f|ImA| 



|A | -> oo . 



So we have 



lim / $/(a;;A)c?A = 2itif{x) 



It is easy to check that 



lim / X)dX = . 



Therefore from the Cauchy's theorem (I3.20p we have 

m - 



^4 n! (2ft 277 1) e _ u /2 u «- n -l/2£2K-2n-l^ 



n=0 



r(2/c-n) 

/■oo 

x / e- tt '/ 2 tt" s - n - 1 / 2 L^- an - 1 («')/(a/)dx' 



+ 



2vr 2 



dA sinh(27rVA) 



r^-ft + zVX) ^ 1/2 W K)i VA(«) 



x / ^ 1/2 ^,Va(«W)^', (3-27) 
with it = 2ft e -( x ~ x o) anc [ similarly for v! . We define the normalized eigenfunctions a£| 



^n(^) 

^a(^) 



n\ (2ft - 2n - 1' 



e -u/2 u ,-n-l/2 L 2 K - 



2n-l, 



r(2ft-n) 
= — sin h 1/2 (2vrv / A) rQ-ft + iVA^ 



u 



u 



(3.28) 

1/2 W_ iVX (u), (3.29) 



with ri = 0, 1, ■■•,[« — 1/2] and A > 0. Then the expansion formula ( 13. 2 7ft leads to the 
completeness relation 



n=0 



(3.30) 



This is the main result of the present paper, which depends crucially on the correct 
normalization (I3.28j) and (I3.29P of the eigenfunctions. The normalization of continuum 
states seems not be discussed in the literature so far. In this section we have derived it 
by the method of contour integral. 



4 Application to the Asian option pricing 

In this section we will apply the spectral decomposition formula (13.271) or (13.301) to the 
problem of Asian option pricing. Let S(t) denote the price of a stock and A(t) = 
| L S(t) dt the time average up to the time t. An Asian put option with the strike 
price K is a contract which gives you the right to sale this stock at the expiration time 
t = T and obtain the payoff K — A(T) if K > A(T), while continue to hold the stock 

4 These eigenfunctions are actually real functions. 
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if K < A(T). Therefore in either case the payoff is P{T) = (K — A(T)) + J§ Obviously 
you should pay some money when signing this option contract at t — 0. Our aim is to 
determine its fair price. 

Usually it is assumed that the stock price S(t) follows the law of geometric Brownian 
motion, i.e. S(t) = Soe crW ^ + ^ r ^ cr2 ^ 2 ^ t , where W(t) is a Brownian motion with constant 
volatility a, and r is the risk-free interest rate. After some rescalings the value of an Asian 
put option can be written as [11] 



-rT 



4So 
a 2 T 



(k-a(r)) + } } (4.1) 



where we have introduced the exponential functional a(t) of the Brownian motion W(t) as 
a{t) = J* e 2 ( w ^ +vt ^dt, together with dimensionless parameters r = <r 2 T/4, v = 2r/a 2 — 1 
and k = tK/Sq. The probability density function f(a,t) of the stochastic process a(t) 
satisfies the following diffusion-type equation 



df d 2 f df 

_£ =2a >_£ + (2(l , + 1 ) 0+1 ).£, a >0. (4.2) 

The heat kernel K(a, a'; t) is the solution of the above equation with the initial condition 
K(a, a'; 0) = 5(a — a'). Then the option value can be written as 



e~ rT 



(^) j™(k-a) + K(a,0-T)da. (4.3) 



Suppose fx satisfies 



_ 2e .0_ (2( „ +1)(l + 1) ^ = 2AA , (4.4) 



which can be rewritten in the form —(pf)' = Xwf by defining p(a) = a y+1 e 1 ^ 2a ^ and 
w(a) = a u ~ 1 e~ 1 /( 2a \ The heat kernel K(a,a';t) has the following spectral decomposition 

K(a,a';t) = w(a) ^ e- 2M f x (a)f x (a') . (4.5) 



Introduce the new variable x and new function ^a^) as suggested in 

x = logo, Va(z) = a^e-^Ma) , (4.6) 
then the equation (I4.4p can be reduced to 

+ \T7 e —e x + — )^ x =X^x- (4.7) 



dx 2 V 16 4 4 



5 We use the definition X + = X- 6(X), where 9(-) is the Heaviside step function. 
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Actually, up to a additional constant z/ 2 /4, this is the static Schrodinger equation with 
the Morse potential. Compared with it standard form ( 13. ip . we have 



1 - v 



1 



2(1- u) 



(4.8) 



Here we restrict to the case v < 1. By use of our results (13 .28j) and (I3.29p in the previous 
section, we can obtain the complete system of eigenfunctions of (14.41) as follows 



fn(a) 
fx(a) 



n\{-V-2n) n j- u -2n ( J_ 



(4.9) 



- sinh 1/2 (2ny/\- v 2 /A^j v(u/2 + iy/X - v 2 /A ) 



xa (.l-u)/2 e l/(4a) w 



(l-v)/2,iy/X-v*/4, \ 2 



1 



(4.10) 



The corresponding eigenvalues are A n = n(—u — n) with n — 0, 1, • • • , [—v/2], together 
with the continuous spectrum A > u 2 /A. Therefore, by redefining p = 2yA — z^ 2 /4, we 
have 



— 

/ — ' rfl — v — n) 



-l/(2a) L 



—v—2n 
n 



n=0 



+ 



2a 



2vr 2 



e -(p 2 +^)r/2^ 2a ^(^l)/2 e -l/(4a) 



X W(l_,,)/2,ip/2 ( ^ 



r 



z/ + ip 



sinh(ir p) p dp . (4-H) 



We have used the asymptotic relation L°(z) ~ - w V z n and ~ e Z / 2 2 K as z — >■ oo. 

Therefore the value of an Asian put option is 



-tT 



o\ P^? (-!)"(-!/ -2n) 



4<Sq \ i \ ^ 

a 2 T/ ] ^ 2n!T(l-z/-n) 

I n=0 



e -2n(-v-n)T Q^y+n+S g -l/(2*)^-i/-2n / _ 

2/c 



+ 



1 Z" 00 

J_ / e -( P 2 +. 2 )r/2 (2A; w+3 e -l/(4fc) 
^ 2 70 



X W_v + 3 ip 

2 '2 \ 2fc 



r 



sinh(7rp)p<ip > . (4.12) 



Here we have used the formula (7.623.7) of [14J to work out the integration over a in (14.31) . 
This is completely equivalent to the result of [13] which is obtained by a different method. 



10 



5 Conclusions 

In this paper we obtain the properly normalized eigenfunctions of the Morse potential 
and explicitly prove the completeness relation by use of the method of contour integrals. 
In this method the discrete spectrum corresponds the pole of the resolvent kernel, while 
the continuum spectrum comes from the contributions of the branch cut. We also apply 
this spectral decomposition to the problem of an Asian option pricing. We expect this 
method can be applied to more subtle spectral problems, e.g. issues about self-adjoint 
extensions, in quantum mechanics. 
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